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Introduction 

Symmetry methods are by now recognized as one of the main tools to attack 
deterministic differential equations (both ODEs and PDEs); see e.g. [3, 5, 9, 
16, 21, 22, 25, 26]. The situation is quite different for what concerns stochastic 
differential equations [1, 13, 14, 17, 20]: here, symmetry considerations are of 
course quite widely used by theoretical physicists (see e.g. [4, 15, 18] in the 
context of KPZ theory), but a rigorous and general theory comparable to the 
one developed for deterministic equation is still lacking, and the attention of 
the community working on symmetry methods should maybe be called to this. 

We would like to quote here early work by Misawa [19], mainly concerned 
with conservation laws for stochastic systems, and on the work by Arnold and 
Imkeller [1, 2] for normal forms of stochatic equations. 

In the following I will report on some work I have done - to a large extent 
in collaboration with N. Rodr' iguez Quintero - on symmetries of stochastic 
(Ito) equations, and how these compare with the symmetries of the associated 
diffusion (Fokker-Planck) equations [6, 8, 23, 24]; further details can be found 
in [12, 10]. Part of this was concerned with providing a suitable definition for 
symmetries of a stochastic equation, so I will also show that the symmetries at 
the basis of the construction of the KPZ equation are actually recovered within 
the frame developed here. 

1 Projectable symmetries 

We assume the reader to be fairly familiar with the concept of symmetry of a 
deterministic differential equation, and the methods to characterize these; so 
we will just recall some formula of use in the following, mainly to fix notation, 
identify the class of transformations we are going to consider, and compare with 
the situation for stochastic equations. 
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Let us first consider the case of systems of first order ODEs in i?", 

- f{x,t) = . (1.1) 

The general form of Lie-point vector fields is in this ease, writing di = d/dx^, 
X = T{x,t)dt +^^{x,t)di. Such a vector field, however, is somehow too general: 
indeed, we would like the change on t not to depend on the value assumed by 
the solution x{t): this is particularly relevant when we think of the ODE as 
describing not the motion of a single particle in R", but that of an assembly 
of particles. We would thus like to keep the special role of time, i.e. require 
T = T{t); in other words, we will consider 

Xo = T{x,t)dt+C{x,t)di . (1.2) 

A similar discussion holds for symmetries of PDEs, first or higher order. 
Here we will deal only with scalar PDEs, i.e. equations for a single dependent 
variable u = u{x,t). In this case we would write in general X = T{x,t,u)dt + 
^^{x,t,u)di + ({){x,t,u)du- On physical terms, it is preferable to consider only 
transformations such that the independent variables are transformed indepen- 
dently of the values assumed by the dependent ones, i.e. such that r and the 
do not depend on u; also, here we are primarily interested in evolution equa- 
tions, and in order to guarantee the time keeps its distinguished role we should 
also ask that r does not depend on x. Thus, in the end we want to consider 
transformations of the form 

X = T{t)^^+e{x,t)^,+ct>{x,t,u)^ . (1.3) 

Definition. The vector fields - and in particular symmetry generators - of the 
form (1.2), (1.3) will he called "fiber-preserving", or simpZy projectable. 

Let us now consider in particular a system of ODEs, which we write in the 
form (1.1). The determining equations for the symmetry generators of this are 
easily determined (these should be compared to the determining equations for 
stochastic differential equations discussed in the next section). We will use, here 
and below, the notation 

:= {P-dj)C-{^'-d,)r . (1.4) 

Lemma 1. The projectable symmetries of (1.9) are given by vector fields Xq 
as above with coefficients satisfying 

5t(e-rr) + {/,a^ = 0. (1.5) 

For given functions f''{x,t) the partial differential equations (1-16) always 
[21, 25] have non-trivial solutions r(t) and Ci^j ^ot example we can fix T{t) or 
^^{x, t) and compute ^*(a;, t) or r(t). Although an infinite number of symmetries 
of (1.16) exist, in general there is no constructive way to find them. 

For further detail on symmetries of (deterministic) ODEs and PDEs, and 
applications, the reader is referred to [3, 5, 9, 16, 21, 22, 25, 26]. 
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2 Stochastic differential equations 



Let us turn to stochastic equations [1, 7, 13, 14, 17, 20], and let us consider an 
Ito equation^ 

dx' = f{x,t)dt + al{x,t)dw^ , (2.1) 

where / and a are smooth functions, a{x,t) is a nonzero matrix and the w'' are 
independent homogeneous standard Wiener processes, so that 

The equation (2.1) should be seen as a map from the vector Wiener process w(t) 
= {'w^{t), w"(i)} to the stochastic process undergone by {x^{t), ...,x"{t)}, 
and its meaning is precisely that of defining the vector stochastic process x(t). 

If wc consider a fmiction y = $(a;), its evolution under the Ito equation (2.1) 
is described (again in terms of a stochastic equation) by the Ito formula: 

"-y a^j";^ 2 aJTaP? "5^ ^ ^ ,n , (2.2) 

= [Pdj<^>' + (l/2)dj„,<i>'{aa^y"'] dt + [{dj(f)')a^''] dw^ . ^ ' 

To the Ito equation (2.1) is associated the corresponding Fokker- Planck equa- 
tion 

dtu = -diifu) + [{aa^ru] , (2.3) 

describing the evolution of the probability measure p = u{x, t) under the stochas- 
tic process described by (2.1). Equations (2.1) and (2.3) contains the same 
statistical information [1, 13, 14, 17, 20] as far as one-particle processes are con- 
cerned (but not if we think the Ito equation as describing a random dynamical 
system, see [1, 10]), provided a satisfies the non-degeneracy condition 

which we will assume throughout this note. 

We stress that we are interested in equation (2.3) as far as it describes 
the time evolution of the probability measure p = u{x, t) under the stochastic 
process (2.1). It is obvious that for this interpretation u{x,t) should be non- 
negative and subject to the condition 



+00 

u{x,t)dxK..dx'' = 1 . (2.4) 



This is relevant in connection with the allowed transformations in the (x, t\ v) 
space: only transformations preserving this normalization do represent symme- 
tries of the Fokker-Planck equation compatible with its probabilistic interpre- 
tation^, and one should expect a correspondence between symmetries of the Ito 

^It is well known that any Ito equation is equivalent to a Stratonovich equation; for the 
Ito equation (2.1), the associated Stratonovich equation is dx^ = 6'(a;,t)dt -|- <T^(a;,t) o dw'^, 
with b' = f — {l/2)[(7{da'^ /dx)]. We will to the Ito formulation in all of this note. 

^In other words, we should set the FP equation in the function space defined by (2.4), 
and thus accept only transformations mapping to itself. 
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equation and this subclass - rather than all - of the symmetries of the Fokker- 
Planck equation. 

Another simple but relevant point which should be stressed is the following: 

different Ito equations which have the same / but different matrices a can 
provide the same term aa^ and thus the same Fokker-Planck equation. 

A simple example is provided e.g. by u orthogonal: in this case we have 
by definition gg^ = I, so that all the Ito equations with the same / and 
any orthogonal matrix a give the same Fokker-Planck equation"^. Similarly, a 
and a = aB, with B any orthogonal matrix, will give the same Fokker-Planck 
equation (and conversely a and a give the same Fokker-Planck with the same 
/ - only if there is an orthogonal matrix B such that the above relation is 
satisfied. 

Thus if we consider a continuous variation of cr, say a + £7, in the Ito 
equation, the associated Fokker-Planck equation remains unchanged provided 
(ct + £7) (a + £7)^ = cra^, which at order e is simply 

(T7^ -t- 70-'^ = . (2.5) 



3 Symmetries of stochastic ODEs 

3.1 Spatial symmetries 

We will at first consider symmetries not acting on the time variable t, i.e. Xq = 
t)di. Thus, we consider a near-identity change of coordinates, passing from 
X to y via 

x' ^y' ^x' +eC{x,t) . (3.1) 

Using the Ito formula, we easily check that at first order in £, this transformation 
maps the Ito equation (2.1) into a new Ito equation 

= f{y,t)dt + al{y,t)dw^ , (3.2) 

where we have explicitly 

f= f + e [dtC + Pd^e - ^'djf + (l/2)(a<7^)^-^a|fee] .3 3^ 

When (3.1) maps (2.1) into itself, i.e. when (3.2) coincides with (2.1) [up 
to terms o{e)], we say that Xq = ^^{x,t)di is a (Lie-point) spatial symmetry 
of (2.1). Thus, Lie-point spatial symmetries of (2.1) are identified by the van- 
ishing of terms 0{e) in (3.3). In the following statement, we use the notations 
introduced above for A and {.,.}. 

•^Thc one-point stochastic processes described by two difEcrent Ito equations having the 
same associated Fokker-Planck equation have the same statistical properties (the probability 
measures evolve in the same way), but are diflferent: the same realization of the Wiener 
process w(t) leads to different sample paths. 
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Proposition 1. The (generators of) Lie-point spatial symmetries of the Ito 
equation (2.1) are given by vector fields Xq = ^^{x,t)dl. with coefficients 
satisfying the determining equations for spatial symmetries of an Ito equation: 

dtC + {/, + {imaa^y^d%C =0 



These arc n + equations for the n + 1 functions r; thus for n > 1 they 
have some solution only in very exceptional cases. This should not be surprising, 
as symmetry is a non-generic property. 

Note that if cr = a{t) does not depend on the spatial variables, the second 
of (3.4) reduces to 'J'^djC = 0) which in turns imply the vanishing of the term 
((Tcr^)^^c??j.^% which in this case can be rewritten as a^dkicr^dj^'^). Thus, (3.4) 
are in this case equivalent to the determining equations for symmetries of the 
deterministic part of the Ito equation [see (1.5) and recall now we are assuming 
T = 0] with the additional condition crj^dj^^ = 0. 



3.2 Spatio-temporal symmetries. 

We could of course consider transformations acting also on t: in this case the t 
transformation would also affect the processes w{t), and some extra care should 
be paid. We will not go over the details (which consist in studying how w{t) 
is changed and include this in the change undergone by a), and just report the 
result obtained in [12], with the same notation as above. 

Proposition 2. The projectable vector field Xq = T{t)dt + ^*(x,t)9J, is a sym- 
metry generator for the Ito equation (2.1) if and only if the coefficients {r,^*} 
satisfy the full determining equations for projectable symmetries of an Ito equa- 
tion: 

dt (e - rf^) + {/, - A^'^d^e = 

{C7k,^y-rdtal-{l/2)aldtr = . ^""-^^ 

Note that the symmetries which are linear in x, i.e. such that ^*(f;x) = 
Mj{t)x^ , are given by the same equations as for symmetries of the deterministic 
part of the Ito equation (2.1), i.e. (1.5), plus the additional condition MjO-j, = 
Mlx^d,ai+Tdtal + {l/2)aldtT. 

A relevant case in applications is the one where /* = /*(x) and (j\ = cr^(t) 
or even al = const = SI; this correspond to an autonomous dynamical system 
subject to a noise which depends only on t or even a constant noise. In this 
case (3.4) can be discussed quite completely. Indeed, with a independent of 
the spatial coordinates, the second of these reads 2cr^9j^* = 2Tdtcrl + cF^dtT. 
As the r.h.s. only depends on i, by differentiating with respect to x™ we get 
the equation al.{d'^£^^/dx^dx"^) = 0: i.e., the (symmetric) matrix of second 
derivatives of Hj^ must be such that aH^ = H^a = 0. Notice that in 
particular, if cr^(t) = ^{i){t)Sl with all A(j) ^ (or however if exists), this 
means that ^ can be at most linear in the x. 
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4 Symmetries of the Fokker-Planck equation. 



We will now discuss the determining equations for projectable symmetry genera- 
tors of the Fokker-Planck equation (2.3) in arbitrary spatial dimensions; general 
symmetries of the Fokker-Planck equations (with some limitations on a) in one 
and two space dimensions have been completely classified [6, 8, 24]. 
It will be convenient to rewrite the Fokker-Planck equation as 

ut + A'^dfjU + B'diU + Cu = ; (4.1) 

the coefiicients A, B, C depend on x and t only and are given explicitly by^ 

A = -(1/2)<T£7^ ; B' = f + 2djA'^ ; C{x, t) = {di ■ f) + djjA^^ . (4.2) 

We consider a projectable vector field of the form (1.3). Computing the 
second prolongation of this and applying it to the equation, we obtain the 
determining equations for symmetries of the Fokker-Planck equation. 

It follows by general results on Lie-point symmetries of linear equations 
[5, 9, 21], or by an elementary explicit computation, that: 

Lemma 2. The projectable symmetries of the Fokker-Planck equation (4-1) o,Te 

given by vector fields in the form, (1.8) with (f> = a(x,t) + f3(x,t)u. 

Again by linearity, we will have "trivial" symmetries = a{x,t)du, with 
a{x, t) an arbitrary solution to the FP equation itself; this is just expressing 
the linear superposition principle [21]. Note that in order to preserve the nor- 
malization (2.4), we should moreover require / a{x,t)dx = 0, which rules out 
nontrivial non-negative solutions. 

We will thus from now on focus on other, i.e. "nontrivial" , symmetries only, 
and take a{x, t) = 0. 

Proposition 3. The nontrivial projectable symmetries of the Fokker-Planck 
equation (4-1) are given by vector fields in the form (1.3) with (j) = (3{x,t)u, 
where r, ^% /3 satisfy the determining equations 

dtirA''') + {£,"'drnA''' - A'^'^dm^'' - A"'''dmC) ^ 

dtirB') - [Ct + B"'drae - e'drnB'] + {A'^dk^ + A'^'dmP) - A"''' d^^C = 
dtirC) +Pt+ A'f'dlP + B'd.P + ^"'draC = . 

(4.3) 

It should be stressed that in the computations leading to (4.3) we have not 
required preservation of the normalization J u dx = 1\ when we require this, it 
turns out that 

Lemma 3. The nontrivial projectable symmetries of the Fokker-Planck 
equation (4-1), in the form (1.3) with 4> = P{x,t)u, preserve the normaliza- 
tion (2.4) if and only if 

P = -div(0 . 

^Recall that our non-degeneracy assumption guarantees A ^ 0. 
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In other words, this lemma guarantees that if we are looking for the (nontriv- 
ial) symmetries of a Fokker-Planck equation compatible with its probabilistic 
interpretation, we do not have to deal with the general form (1.3) of the sym- 
metry vector field, but we can use instead the ansatz 

X = T{t)dt + e{x,t)di - [diy{Ou]d^ ; (4.4) 

this substantially simplifies the analysis of symmetries. 

5 Symmetries of the Ito versus symmetries of 
the associated FP equation. 

We are specially interested in discussing how the symmetries of the partial 
differential equation (4.1) and those of the symmetries of the system of stochastic 
ODEs (2.1) are related. 

In order to do this, we should express the coefficients A, B, C in terms of 

/ and a; indeed A = -(l/2)crcr^, = .f + 2dkA''' , and C = 9,.f + df^A''' . 
Using the latter two of these and with some manipulations, see [12], (4.3) reads 

dtirA''') + {C^'dmA''' - A'"'d.mi'' - A''"'drr,e) = 

dtiC - rf ) + {/, - A"^^dl,e] - 2 [A^'dkP + A^^^dl^e] = (5.1) 

dt + fdi - A^f^dl] [(3 + dmD = . 

It is convenient to introduce the shorthand notations 

r| = aj^d^e - rdm<T^ - Tdta^ - il/2)a'jdtT , 

A' = - [dtiC - rf) + {/, - ^'"'a^fef] • ^ ' 

With these, (5.1) is rewritten as 

crjrjjj^ + a^ris^^ = 

A' + 2[A'''dkP+ (5-3) 
[dt + fdi - A^'^dl] [(3 + 9„r ] = . 

It should be stressed that with the notation just introduced, the determining 
equations (3.5) for symmetries of the Ito equation (2.1) read simply 

A' = ; rj = . (5.4) 

With this notation, the comparison of symmetries for the Ito and the asso- 
ciated Fokker-Planck equation is quite a simple matter. 

5.1 From the Ito to the FP symmetries 

We will first investigate if symmetries of an Ito equation result in symmetries 
of the associated Fokker-Planck equation. 
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Wc first note that for F = 0, the first of (5.3) is satisfied. As for the second 
of (5.3), when A' = this reduces to 



(5.5) 



Obviously [3 is not present in symmetries of the Ito equation, so we can 
choose it as to satisfy the third of (5.3) and (5.5): for this it suffices to choose 



We know from lemma 2 that actually wc have to take co = 0, and that this 
is actually the only choice of (3 providing us with a normalization-preserving 
symmetry of the FP. Our findings are summarized as follows: 

Proposition 4. Let Xq = rdt + C^i be a symmetry of the Ito equation 

(2.1). Then Aq extends to a unique normalization preserving symmetry X := 
Xq — [{div^)u]du of the associated Fokker-Planck equation. 

5.2 Prom the FP to the Ito symmetries 

Let us now consider the converse question, i.e. if and when a symmetry of the 
FP equation associated to an Ito equation can be projected to a symmetry of 
the Ito equation itself. 

As the probabilistic interpretation of the FP requires the normalization con- 
dition (2.4), we expect only symmetries identified by lemma 2 can provide sym- 
metries of related Ito equations. 

Moreover, as recalled above, there are transformations which map an Ito 
equation into a different one with the same statistical properties and thus the 
same associated FP equation: these would be symmetries of the FP but not of 
the Ito equation. Thus, not all the symmetries of the FP should be expected to 
produce symmetries of the Ito equation. 

By lemma 2, we assume /? = — div(^); this guarantees that the third of (5.3) 
holds, and also that the second of (5.3) reduces to A* = 0. Hence we only have 
to discuss the relation between the first of (5.3) and the second of (5.4). These 
two equations can be rewritten, respectively, as 



The second of these imply the first, but the converse is not true: hence there is 

not a complete equivalence. 

This just corresponds to a given FP equation corresponding to different Ito 
equations: the transformations with F ^ but aV^ + T<j^ = will be precisely 

those which map an Ito equation Ea into a different Ito equation E\ which has 
the same Fokker-Planck associated equation, see (2.5). 

Proposition 5. Let X = T{t)dt + S,^{x,t)di — [(div^)u](9„ be a symmetry 
of the Fokker-Planck equation associated to the Ito equation (2.1). Then X 
maps the Ito equation into a (generally, different) Ito equation with the same 



P = -dmr + Co = -div(0+C0 . 



(5.6) 



aV'^ + Ta^ = and F = . 



(5.7) 
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statistical properties; X is a symmetry of the Ito equation (2.1) if and only ifV 
defined in (5.4) satisfies = 0. 

We may note that from the above system (5.4), restricting to the case where 
f{t,x) = f{x) and (j{t,x) = <t{x), we recover the results of Cicogna and Vitali 
[6] for the one- dimensional setting. From (5.2) it is also possible to recover the 
results of Shtelen and Stogny [24] for the two-dimensional Kramers equation, as 
well as the results of Finkel [8] for certain two-dimensional FP equations. 



5.3 Examples 

Example 1. As the first example, we consider the case f{t,x) = 0, a{t,x) = 
(To = const ^ 0, i.e. the equation 

dx = (To dw(t) (5-8) 

which represents a free particle subject to constant noise. The corresponding 
Fokker-Planck equation is simply the heat equation ut = {aQ/2)uxx- The sym- 
metries of the heat equation (other than the trivial ones Va = a{x,t)du with 
Oit = oixx) are well known to be [21, 25] 

vi = dt ; V2 ^ dx ; vs = udu ; W4 = (y^tdx - a^xudu ; ,r. ^-^ 

V5 = 2tdt + xdx ; W6 = t'dt + xtdx ' {l/2){t + x^/al)udu ^ ' 

Of these, v\, and (which do not act on ti) are also symmetries of the Ito 
equation (5.8), as is easily checked using (3.4). Notice that (4.4) is satisfied for 
these, and is not satisfied for W3,W4 and wg. The vector fields Vi,V2 and do 
actually span the symmetry algebra of (5.8). 

Example 2. As an example in two space dimensions [with coordinates 

{x\,X2) = {x,y)], we choose 



dx = ydt 

dy = -k'^ydt + V2k'^dw(t) 



with a positive constant. 

The corresponding Fokker-Planck equation is the Kramers equation 

Ut = k^Uyy — yux + k^yuy + k'^u ; (5.10) 

the symmetries of this were studied in [24] and, apart from the trivial ones Va, 
are 

vi =dt ; V2 = dx ; w-s = e"''^* [k~^dx - dy] ; V4, = udu ; 

V5 = tdx + dy- (1/2) (y + Px)udu ; vq = e^^^ [k~^dx + dy - yudu] 

Here vi, V2 and 113 satisfy (4.4), while for V4, W5 and vq this is violated. According 
to our definition, the symmetries of the equations (5.9) are again vi, V2 and V3, 
as easily checked. 
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Example 3. Wc consider now an example where the correspondence be- 
tween normahzation-preserving symmetries of the FP equation and symmetries 
of the Ito equation is not complete, i.e. the two-dimensional Ito system (with 
zero drift) 

dx^ = cos{t)dw^ — sm{t)dw^ , 
dx^ = sm{t)d'W^ + cos(i)rfw^ ; \ ■ ) 

the corresponding Fokker-Planck equation is now just the two dimensional heat 
equation Ui = (1/2) A?!. 

It is now immediate to check that the vector field Xq = dt is a symmetry 
of the FP equation, but not a symmetry of (5.11). Obviously, the case of any 
orthogonal a with dtcrl ^ will be exactly the same. 

Example 4. Consider n uncoupled equations for equal "Langevin harmonic 
oscillators" subject to independent stochastic noises [17, 23]; this system is 
described by the Ito system 

dx' = -x" dt + y/2sidw\ i=l,2,...,n (5.12) 

(no sum on i) where we assume all the Sj are strictly positive; the corresponding 
Fokker-Planck equation is 

n 

dtu = ^ [sidfiU + x'diU + u] . 

i=l 

In this case we get the symmetries (i = 1, 2, n) 

n 

vi=dt , V2 = e"^*[5t -^x'di + nudu] , % = e"*5j . 



Exemiple 5. Finally, we consider now the n-dimensional nonlinear case 

dx' = -(1 - \\\x\\^)x' dt + dw' , (5.13) 

where ||a;|| is the norm of the vector x and A ^ 0. By inserting /* = — .x*(l — 
A||a;||^) and a] — 5ij in (5.4) we obtain, from the second of these, ^' = h^{t) + 
(Tt/2)a;'. Inserting this into the first of (5.4) and isolating the coefficients of 
different powers of x, we get that the only symmetry is given by v = dt- 

Note that we got no rotation symmetry; this is because within the class 
of transformations we are considering we can rotate the vector x, but not the 
vector Wiener process w{t); transformations allowing to rotate 'w{t) as well will 
be considered in the next section. 
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6 W-symmetries 



We will consider symmetries involving not only the spatial and time variables 
(x.t), but also the vector Wiener processes w{f,) entering in the n-dimensional 
Ito equation (2.1). We will specifically consider infinitesimal transformations of 
the form 

=x^ + e^\x,t) 
s =t + eT{t) (6.1) 

—^ = + eiJL^{w, t) . 

We also call symmetry generators of this form, " W-symmetries^^ . 

Note that this is not the most general possible form of a transformation for 
the variables involved in (1.1); some words on this restriction are in order here. 

The restriction on r and ^ is the same as considered above. Moreover, we 
only allow "internal" transformations of w, i.e. they cannot depend on x{t). As 
we think of the stochastic process w{t) to be independent of the evolution of the 
x{t), we like its transformation not to depend on the latter. Finally, we have 
allowed the transformation on the spatial coordinates x{t) to depend on x and t, 
but not on the w; this means that we do not want to consider transformations of 
the spatial coordinates which depend on the realization of the stochastic process 
w{t). This again is somewhat a natural requirement in physical terms^. 

In order to compute how (6.1) acts on the Ito equation (2.1), we use the 
Ito formula (2.2). Moreover, due to the form of (6.1), we can adopt a "two- 
steps procedure", see [10]: the transformations of the {x,t) variables and of the 
stochastic process w{t) do not interfere with each other. It turns out that some 
strong limitation on the functions /u' arise. 

Lemma 4. The infinitesimal transformation w ^ z = w + e/i(w,t) maps 
the vector Wiener process w{t) into a vector Wiener process z{t) if and only 
if H = Bw, with B a real antisymmetric matrix. Equivalently, if and only if 
z = Mw with M orthogonal. 

Note that as B is antisymmetric matrix, it vanishes in dimension one: for 
one-dimensional Ito equations we have no new symmetries with respect to those 
discussed above. 

Proposition 6. Under the infinitesimal transformation (6.1), with B a real 

antisymmetric matrix, the Ito equation (2.1) is changed, into a (generally, dif- 
ferent) Ito equation. The Ito equation (2.1) is invariant under (6.1) if and only 
ifTj^j^ satisfy the determining equations for W-symmetries: 

dtc + {f,o' - W) + A^'^d^.e = 



^It should be however stressed that in this way we are also discarding the transformations 
needed to obtain normal forms of stochastic differential equations, see [1, 2] and references 
therein. 
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If we just consider a transformation w ^ z = Mw, with M orthogonal as 
required by lemma 4, the Fokker-Planck equations associated to the Ito equation 
and to the transformed one will be the same. Thus all the discussion conducted 

in section 5 about the relations bctwcn symmetries of an Ito equation and of 
the associated FP equation also applies to transformations of the form (6.1)/' 

Example 1. Let us consider again the two-dimensional Ito equation (5.9) 
If we look at the second set of determining equations and single out the one for 
i = 2 and fc = 1, we have immediately that 5 = 0, i.e. in this case we have no 
new symmetry by allowing transformation of w. 

Example 2. Let us consider again the two-dimensional Ito system (5.11) 
We will write = a, S,^ — (3, and B12 = —B21 = b. The most general solution 
to the determining equations turns out to be 

a = [Tt/2-sm{2t)T]x + [b + cos{2t)T]y + aoit) 
13 = -[h + T]x+[Tt/2\v + Po{t) ; 

By setting r = cto = /3o = 0, we obtain a new symmetry, corresponding to 
a = by and /3 = —bx with b arbitrary, i.e. given by 

X = ( — - x—\ + (w^— - w^—\ 
\ dx ^ dy ) \ dw^ ^ dw'^ J 

Note this is nothing else than a simultaneous (and identical) rotation in the 
{x,y) and in the (wi,«;2) planes. 

Example 3. We will now consider again the case of n uncoupled equal 
"Langevin harmonic oscillators", see (5.12) above. 

If we discard symmetries with B = (i.e. those obtained above), we are 
still left with a nontrivial possibility, i.e. B an arbitrary (real, antisymmetric) 
constant matrix and = C^^x^, where C\ = [si/skjB]. (no sum on i,k here). 

The meaning of this result is obvious: we can act on this system by an 
arbitrary SO{n) rotation in the w space and a related rotation in the x space. 

If we assume si = S2 = ■•• = Sn = s, we have indeed C = B and the rotations 
in w and x spaces do just coincide. For general si, ...,Sn, the (si/sk) factors 
relating the C and B matrices are also easily understood: wc could rescale each 
of the by a factor y/2si, = ^2siy^ (no sum on i), arriving at the manifestly 
rotationally invariant n-dimensional Ito equation dy* = —y'^dt + dw^. 

Example 4. Let us consider again (5.13). We easily check that it is rota- 
tionally invariant with a W-symmetry: take r = and ^* = B^x'^, so that 

{^jy = {Bx,il-X\\xf)x} = -2\x'{x^Bix^) = 0, 

the last equality following from B = —B^. That is, (5.13) is indeed symmetric 
under simultaneous identical rotations of the x and w vectors. 

^Needless to say, for more general transformations, for which the transformations of the 
{x,t) variables and of the vector Wiener process w{t) do interfere with each other, this would 
not be the case. 
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7 Discrete symmetries 



We will now briefly consider, for the sake of completeness, discrete symmetries 
of stochastic differential equations. Similarly to what happens in the determin- 
istic case (sec e.g. [11]), the resulting determining equations are in general too 
difficult to be attacked except for very simple classes of transformations (e.g. 
reflections, or however linear ones). On the other side, they can be used to check 
if a given discrete transformation is a symmetry of a given stochastic differential 
equation; and/or to determine the stochastic differential equations which admit 
a given discrete transformation as a symmetry. 

We consider again an Ito equation of the form (2.1); here we will not consider 
transformations acting on the t coordinate, and will thus limit to consider the 
change of coordinates in the x and w spaces given by 

(7.1) 

As required by lemma 4, we take R E 0{n). We write as usual A = — (l/2)crcr^. 
With the Ito formula we have at once that y obeys the Ito equation 



y' = <f>'{x,t) , w'^{t)^R%{t)zP{t) 



d4> 



if 



92 



dxidx^ 



+ 



dt 



dt + 



dxo 



.R\ 



dz" 



(7.2) 



Proposition. The transformation (6.1) is a symmetry of the Ito equation (2.1) 
if and only if (j), R satisfy the determining equations for discrete symmetries of 

an Ito equation: 



^-^f^{x,t) + A^fc(x,i)5£^ + 
"liBll„3^i^^^t)R\ = a\{^{x,t),t) . 



r {<^{x,t),t) 



(7.3) 



dx3 



Example 1. The simplest case of discrete transformation is provided by 

(j)^{x,t) = —x^, R = ±1. In these cases, (7.3) reduce to 

fix, t) = -fi-x, t) ; a{x, t) = T<y{-x, t) . (7.4) 



Example 2. For the n independent "Langevin oscillators", see (5.12), it is 
immediate to check that (7.4) is satisfied (for R = —I). More generally, (7.3) 
read now (no sum over i, sum over j when repeated) 

Example 4. For equations with a^j.{x,t) = So^\ (with sq a real constant), 
the second of (7.3) always reduces to {d4>'^/dx^) = (R^^Yj. But this implies 
must be a linear function of the x, (j)'^{x,t) = L^j{t)x^ with L an orthogonal 
matrix [L = R~^). Thus in this case we are always reduced to the simple 
equation 

L'jf{x,t) = f\Lx,t) [LGO(n)] . 
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8 Symmetries of the discrete KPZ equation 



The KPZ (Kardar-Parisi-Zhang) equation is written, in physicists' notation, as 

^ = „Ah + ^ {Vhf + n . (8.1) 

Here h = h{y,t) is a scalar variable, representing the height of a growing (non- 
overhanging) interface at point y € i?" at time t, the operators V and A rep- 
resents the gradient and laplacian with respect to the y variables, v > Q and 
A > are real numbers, and rj ~ 77(7/, t) is a random field with 5 spatial correla- 
tion function and the statistical properties of a white noise at each given point 
y. When A = the KPZ reduce to its linear part, which is also known as the 
Edwards- Wilkinsons equation. 

The KPZ is a universal equation describing, under general hypotheses and 
up to an asymptotic procedure, the time evolution of a growing interface; see 
[4, 15, 18] for details. 

We will limit ourselves to the one-dimensional case n = 1; moreover, as our 
discussion deals only with stochastic ODEs (as opposed to stochastic PDEs), 
we consider a discretization of the KPZ, obtained by considering a imiform 
chain of intersite distance, {5y); in other words, we consider a discrete set of 
points yk = k{6y). Thus h{y,t) is now replaced by a vector variable x{t) with 
components x'^{t) (i E Z), where .T*(t) w h{yi,t). 

The gradient and Laplacian terms in the KPZ will then be replaced by 

m{y,,t) « ; {^h){yut) « ^ • (8.2) 

As for r]{y,t), its discretization will be a standard vector Wiener process w{t) 
with component w'(t) at the point yi. We also write a = i'[2{Sy)]~^ and /3 = 

In this way, we replace (8.1) by the Ito equation 

da;^ = (a[x'+'^ -2x' +x'-'^] + p [x'+'^ - x'-'^f^ dt + dw\t) . (8.3) 
It will be convenient to rewrite this in a more abstract notation as 

dx* = [M\x'' + T\jkX^x''] dt + dw\t) (8.4) 
(note we have a = I); by comparing (8.3) and (8.4) we have 



M\ = a 



(8.5) 



Note that F^j, = F^^; it will be convenient to also write {G(j))\. = F*^,. 

Physically an equation describing a growing interface is required to possess 
the following symmetries [4, 15, 18]: (1) time translation; (2) translation in the 
h direction; (3) space translations; (4) spatial inversion across any plane and 
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(for n > 1) rotations around an axis perpendicular to the x plane. Moreover, 
at the linear level we should also require invariance under: (5) inversion in the 
h directionJ 

It is well known that the KPZ is the simplest nonlinear (stochastic PDE) 
equation possessing these symmetries, and actually it was determined essentially 
on the basis of this consideration. 

We will now check that with our definition of symmetry for a stochastic 
differential equation, the discretized KPZ (8.3) does indeed admit these sym- 
metries. Note that in the discretized case, (1) and (2) remain continuous sym- 
metries, while (3), (4) and (5) are discrete ones. In particular, (3) is replaced 
by the shift of one site along the chain, and (4) reduces to spatial inversion. 

8.1 Continuos symmetries 

We start by analyzing continuos symmetries. The second of the determining 
equations (6.2) provides, for a = I, 

which we also write as ^' = A\x'' + a*. Plugging this into the first of (6.2), and 

separating terms homogeneous of different degree in the x, we get (using also 
symmetry of the F and our shorthand notation) 

ai - Mja^ = 

i [T„5j - (Bt)]] - TtM] + [A, M];. - 2af G(,); = (8.6) 
TtG(j)'j,- A^G(j)^ + G(j)^A^ = 

We will not try to determine the most general solution to these, but just 
check that (l)-(5) are indeed solutions to these equations, i.e. that our defini- 
tions correctly identify them as symmetries of the (discretized) KPZ equation. 

As for (1), in this case r = 1, while = 0, which of course imphes a and F 
are also zero, and B = as well. It is immediate to check that (8.6) are indeed 
satisfied in this case. 

For the transformation (2) we have r = 0, i? = 0, = 1; the latter means 
A = 0, a' = 1. The first and third of (8.6) are immediately seen to hold, while 
the second reduces to Ylj ^ with the explicit expression for F this is also 
easily seen to hold, as it reads 

+00 

E [ci+Ci-ci-ei; • 

j=-oo 

8.2 Discrete symmetries 

In the case of discrete symmetries (3), (4), (5), we deal with linear transfor- 
mations, i.e. a; — > y = Fx, w ^ z = Rw. The relevant equations (7.3) are 

Checking that this is the case for the KPZ equation is immediate. 
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then 

F;f = f{Fx) ; FX{ir')l = 4 . (8.7) 

As a = S, the second of these actually reads R = F. Writing / as in (8.4) and 
separating linear and quadratic terms in the x, the first of (8.7) reads 

3 km pq k m 

We thus have to check these are satisfied, with M and A provided by (8.5), for 
(3),(4),(5). 

In the case of translations, (3), we have = <5*_|_i. Hence 

{FM)i = di-^\ - 2Si-' + Si-_\ , {MF)i = - 2SU, + Si , 
and the first of (8.8) is satisfied. Similarly, 

i^r™ = + si-_\ + <5};1<5^ - 6p, - > 

and the second of (8.8) is satisfied as well. 

Let us come to (4), i.e. spatial inversion with respect to any site m: this is 
described by F^ = (no sum on m in this and the following formulas). In 

this case we have 

{FM)] = <5j^_, - 25;- + 5;-+, ; {MF)) = S^^I^f' - 2<5j^ + . 

This shows that the first of (8.8) holds. As for the second, it also holds as it 

reads 

"m.-j+l"m-k+l ~^ "m-j~l"m-k-l "m-j + l'^m-k-1 "m-j-l"m-fe+l ~ 

ci— m+1 ci— m+1 , ct—m—lci—m—l ci — m+lci—m—l ci— m— 1 ci — m+1 

^m—j m — k ' m—j m — k m—j m — k m—j m — k 

Finally, let us consider (5): this corresponds to F = —I, so obviously 
[F, M] = 0. The second of (8.8) reads in this case F = 0, which only holds 
for A = 0, i.e. in the linear (Edwards- Wilkinson) approximation. 

We conclude that our definition correctly captures the symmetry properties 
of teh KPZ equation. 
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